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We apply the muffin-tin effective medium approximation to calculate the temperature dependence 
of the resistivity and of the thermopower of amorphous and liquid metals. The results show unam- 
biguously that a large resistivity is accompanied by a negative temperature coefficient, in agreement 
with the experimental situation. This behavior is shown to result from a pseudo-gap which opens 
in the 1-particle spectrum due to strong scattering at the quasi zone boundary. In turn the ther- 
mopower is found to have non-trivial density and temperature dependences. 
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I. INTRODUCTION 



The influence of a disordered potential on the motion 
of electrons remains an interesting problem despite of 
decades of work. The aim of solving this problem is 
to elucidate the properties, if not the existence, of a 
whole range of materials, such as doped semi-conductors, 
doped Mott insulators, amorphous and liquid metals, to 
quote a few. The classes of materials are not only distin- 
guished by the state of their parent compounds, but also 
by the type of disorder and the influence of the electron- 
electron interaction. Even if the latter does certainly 
play a key role for the doped Mott insulators, it is com- 
monly accepted that its role is much reduced in most 
amorphous and liquid metal systems. Types of disor- 
der can be, among others, "diagonal" as in conventional 
modeling of doped semi-conductors, which can leadr-to 
the localization transitional (for a recent review seeB), 
or topological, or both. The strength of the disorder is 
an equally important characterization of a physical sys- 
tem. In amorphous, and liquid, metals the disorder is by 
far not extreme. Instead, neutron scattering data reveal 
distinct short ranged order (for a review, seeS) and the 
ionic structure factor displays a series of peaks. Never- 
theless, those metals exhibit a whole range of anomalous 
properties; especially transport properties exhibit pecu- 
liar behavior, and even the electronic structure poses a 
challenge to the theoretician. How to explain the mini- 
mum in the Density of States in noble metal-polyvalent 
metal glassesQ orreven the true gap in the spectrum of 
liquid Bi, Tl, Hga. What is the influence of the elec- 
tronic structure on the transport? On the experimental 
side, detailed analysis of noble metal-polyvalent simple 
element glassesQ have revealed a strong influence of both 
the ionic structure on the electronic one, and of the elec- 
tronic structure on the ionic one. Namely when K p (the 
location of the first peak of the structure factor) is about 
2kp (the diameter of the Fermi-sphere), a minimum de- 



velops in the electronic density of states, which is best 
seen in these particular alloys. In turn it causes the 
Friedel oscillations of the effective electron-electron inter- 
action to match the ionic potential (at least as far as the 
first few minima are concerned), thus tending to stabi- 
lize the amorphous structures in a..manner similar to the 
cij^stalline Hume-Rothery phasesffl. Recently Kroha et 
aS3 obtained an instability of the Fermi see towards den- 
sity oscillations which causes a phase shift of the Friedel 
oscillations, as well as an enhancement of their ampli- 
tude. This instability occurs when the transport time 
is substantially smaller than the quasi-particle life-time 
and is consistent with experiment. Experimentally this 
behavior is best seen when the valence of the glass Z 
is close to 1.8 el/Atom. In these glasses transport co- 
efficients have been systematically investigated too. As 
a function of valence, the electrical resistivity is maxi- 
mum for Z ~ 1.8. There it decreases with temperature, 
while the thermopower grows with the temperature at 
low temperature. The correlation minimum in the den- 
sity of states, negative temperature coefficient of the re- 
sistivity and positive thermopower has not yet found a 
definite theoretical explanation. 

To tackle this problem we apply the MufEn^Tin Ef- 
fective Medium Approximation (MT-EMA)Eftd, to an 
ensemble of atoms, distributed according to the hard 
sphere structure factor. Having the electrons moving in 
the muffin-tin potential created by the atoms is an ap- 
pealing description of those all-in-all rather good metals. 
The MT-EMA is a self-consistent scheme which is re- 
summing the disorder averaged multiple scattering series 
and addresses the 1-particle and 2-particles quantities on 
an equal footing, so as to fulfill the Ward identities. It 
precisely allows for a detailed discussion of the interplay 
of multiple scattering and short ranged order effects. A 
first study existsEj. Here we extend it to a complete dis- 
cussion of the temperature dependence of the electrical 
resistivity and of the thermopower. 
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II. SUMMARY OF THE FORMULATION 



In this work we investigate the effect of short ranged 
order (SRO) in an otherwise topologically disordered sys- 
tem, on the electronic structure and transport properties 
of the medium. Such a question is relevant to metallic 
glasses and liquid metals, where the periodic arrange- 
ment of the ions is replaced by a structure characterized 
by SRO only. In this case a series of properties of the met- 
als are anomalous. Most striking is the behavior of the 
electrical conductivity as a function of the tempgrature 
A compilation of experimental data by MooijEElO leads 
to the " Mooij rule " and separates the alloys into 2 classes 
of behaviors: first if the low-temperature resistivity (po) 
exceeds a certain threshold value p c (about 150 pflcm), 
then increasing the temperature leads to a decrease of 
the resistivity, and second if po < p c the more common 
behavior of seeing the resistivity increasing upon an in- 
crease of the temperature is restored. In the meantime, 
a series of exceptions to the second class of behavior has 
been discovered^, and only the first behavior holds for 
good. Quite remarkably p c is not particularly large, not 
even 2 orders of magnitude larger than in good metals 
such as Cu at room temperature. It nevertheless corre- 
sponds to a strong reduction of the mean free path down 
to several atomic spacings. In this regime weak scattering 
approaches are thus expected to be inappropriate. 

Another interesting quantity is the thermopower S. 
While at low temperature a simple estimate using Boltz- 
mann equation yields S/T = —k 2, n 2 l\e\Ep, experimen- 
tal evidence points towards a very different situation. 
Not only the magnitude, but even the sign of S/T are 
not universal. They strongly depend on the density of 
charge carriers. The Hall conductivity and the magne- 
toresistance are anomalous, too, but their investigation 
is left for future work. 

Given the short mean free path, the elucidation of 
these anomalies is not expected to follow from a sim- 
ple weak scattering approach, such as Ziman ihrmula for 
the resistivityEa Indeed, as noted by Jacklej23, the am- 
plitude of S/T can be too large to be obtained out of 
the Ziman formula, even though the general features can 
be understood in this framework. Moreover the exper- 
imental finding of a structure induced minimum in the 
density of states of noble metal- polyvalent metal glasses 
by Haussler et al forces us to abandon the nearly free 
electron model. Instead we resort to the MT-EMA. In 
this framework free electrons are scattered by the po- 
tential of the ions, and the self-energy (which allows for 
interpreting the results in terms of an effective dispersion 
relation) is determined by a self-consistent scheme which 
takes SRO effects into account, by resumming approxi- 
mately the disorder averaged multiple scattering series. 
This scheme becomes exact in the limit of large coordina- 
tion numbee, provided the Kirkwood superposition prin- 
ciple holdsEil. In real systems the coordination is fairly 
large (about 11), and allows us to conclude that we have 



quite a reliable approximation. In turn transport coeffi- 
cients are determined self-consistently, such as to fulfill 
Ward identities. 

We restrict ourselves to the case of s-scattering only. In 
sum, we first determine the (disorder averaged) Green's 
Function: 

<G k (E) > = Go(fc) + G ok (E)Q k (k,k)G ak (E) (1) 

where the off-shell part of the (disorder averaged) T- 
matrix is related to the on-shell part at k 2 = E by: 



Qk(p,p') = n 



t(p,p')r - t{p, «)*(«, j/) 



Qk{K, k) 



(2) 



the on-shell part of the T-matrix is obtained in terms of 
its diagonal part td and an effective propagator G k (k) as: 



nt d (K, k) 



1 - nt d (n, n)Gk(n) 
The latter 2 quantities are obtained as the solutions of 



(3) 
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a (k, q)G q {n)Q q {K, K)G q {n) 



t d (n, K) = t(K, K) 

/k 2 dk ~ 
-^ r G k {K)Q k {K,K)B k {K)t d {n,K) (4) 

where we have introduced t(jp,p') as the ion off-shell 
t-matrix, B k (n) the effective propagator in the Quasi- 
Crystalline Approximation 



B k ( K ) = B fc (/c) 



q 2 dq 
2it 2 n 



a (k 1 q)B q (K) (5) 



with B k (n) = 1/(k 2 — k 2 ) being the free electron prop- 
agator. The information about the structure enters via 
the pair distribution function h(k) (which is related to 
the structure factor by S(k) = 1 + h(k)) and the density 
of the ions n. Expending the pair distribution function 
according to: 



h(\k - q|) = 4TT^Y L (k)ai(k,q)Y L (q) 



(6) 



defines the quantity ao(k, q). We obtain the density of 
states in the usual way: 



(7) 



mr J {2-K) i 
and the electronic density by 

n e = [ dEN(E)f F (E-p) . (8) 



Next we solve the equations for the transport. The vertex 
corrections are determined by solving: 
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out of which we construct: 
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which allows for obtaining the E— dependent conductiv- 
ity 

a{E ) = A(^ )2 f ^f Gok{E+)Gok{E -)A k (E) 



- ^Re(Gl(E+) 



d^ k {E^ 
dk 



h ,eh. 



3n m J \ 7r 

out of which the conductivity results into: 
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We shall compare the EMA conductivities to the ones 
following from Ziman formula 
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In linear response the thermopower S results into: 
For more details we refer toBil. 



(13) 



(14) 



III. RESULTS 

We first describe how the atomic structure is taken into 
account. Here we resort to the Percus-Yevick solution of 
the hard sphere liquid to model the ionic structure. The 
resulting structure factor is displayed in fig. [l], for pack- 
ing fractions r\ ranging from rj = .38 to 7/ = .42, which are 
reasonable values in order to account for at least liquid 
simple metals. 




FIG. 1. The hard sphere structure factor for different 
packing fractions. 

Even though merely a model structure factor, it is ex- 
pected to capture the main features of the structure, and 
has beaea favorably compared to the one of liquid simple 
metalsc3. It also has the advantage of substantially re- 
ducing the numerical effort since h(k) has an analytical 
expression. Here the height of the first peak of the struc- 
ture factor can be varied by tuning the packing fraction, 
in order to mimic a change in temperature. In order to 
model the dependence of 77 on T, we have reported the 
height of the first peak of the structure factop-|S'(A:p) as 
a function of temperature for various systemscJin fig. |^. 
Applying a best fit procedure reveals that Srikp) obeys 
a law: 



Sr(k p ) = 



1 



aT + b 



(15) 



sticking to liquid Rb for the sake of definiteness yields 
a=5.09 lO^K- 1 and b = 0.23378. We can then model 
the temperature dependence of rj(T), listed in table 1. 

TABLE I. Temperature dependence of the packing fraction 



S(k p ) 


2.208 


2.129 


2.049 


1.970 


1.901 


V 


.42 


.41 


.40 


.39 


.38 


T[K] 


430 


463 


499 


537 


573 
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We use the square well potential to model the screened 
Coulomb potential of the ions. In atomic units it is de- 
scribed by the parameter V z = Va? , where V is the depth 
of the potential and 2a its range. Throughout this work 



we use V z — 1.4 and the unit h — 2m 
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FIG. 2. 
factor S(k p 



The amplitude of the first peak of the structure 
) as function of temperature. [29] 1: liq. Al (6), 



2: liq. CuSn (1), 3: liq. GeTe (2), 4: liq. AlGe (3), 5: liq. 
Ga (4), 6: liq. Sn (4), 7: liq. Cd (4), 8: liq. Zn (4), 9: liq. 
Zn (5), 10: liq. Rb (7), 11: liq. Cs (8), 12: liq. Ga (9), 13: 
liq. Rb (10), 14: amor. MgZn (11), 15: amor. CaZn (11), 16: 
amor. MgCu (11), 17: liq. Zn (12), 18: liq. Rb [23]. 



A. Electronic structure 

We now proceed to solve the set of equations (|4|) , de- 
termining the electronic structure. To some extent these 
equations have been solved by Fresard et aO. Here we 
solve them in order to first determine the parameter 
range where the theory is analytic. Indeed the break- 
down of analyticity of the MT-EMA has been observed, 
for the delta- function potential by Singh and Rothcll 
possibly yielding negative spectral functions especially 
for large packing fractions, and no proof, of analyticity 
of the approximation has been founcO'Ell It turns out 
that the analyticity breaks down for low Fermi energies, 
when rj > 0.44, namely a value which is above those we 
are dealing with in this work. For r\ < 0.42, the the- 
ory is analytic for all values of the potential and of the 
Fermi energy. The main result concerns the density of 
states and is displayed in fig. |[ 2 features appear clearly. 

First a deep minimum, and second quite a sizeable 
peak just below it. The minimum is becoming deeper 
and deeper as r\ is growing and occurs at E m , with 
E m ~ (kp/2) 2 with little dependence on the potential 
strength. It is thus fair to say that it results from quasi 
Bragg scattering on a quasi zone boundary. It replaces 



the energy gap which opens when the lattice periodicity 
is restored. As a result breaking the periodicity of the 
lattice does not invalidate the concept of zone boundary 
as a whole: even though the latter cannot be defined as 
the inverse of the lattice spacing but is rather given in 
an effective way by about k p /2, it still locates the region 
in k-space where a pseudo-gap opens due to strong back- 
scattering. At E m , and below it, the effective dispersion 
flattens which leads to the maximum in the density of 
states. Such a maximum leads to a decrease in the total 
energy making it very likely to find most real liquid and 
amorphous metals with a Fermi energy close to r E' m , in 
agreement with the argument of Nagel and Taucn3. 
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FIG. 3. The electronic density of states, for the potential 
strength V z = 1.4 and for different packing fractions. The 
dashed curve corresponds to the calculated free electron den- 
sity of states. Inset: expanded view around E m . 
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FIG. 4. The density of electrons, for the potential strength 
V z = 1.4 and for different packing fractions. Inset: expanded 
view around E m . 
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The dependence of the charge carrier density on the 
Fermi energy is displayed in fig. ^ , and it appears that 
E m corresponds to n m = 1.78. This value is not uni- 
versal, and depends on both V z (clearly) and rj (weakly). 
Typically it will be shifted down to lower densities by 
reducing V z , or increasing rj. Note that the lower band 
edge is shifted down from due to the averaged attractive 
potential. More details concerning the spectral function, 
the-sefl-energy and the effective dispersion can be found 
iiOL 



B. Electrical Resistivity 

We now proceed to the solution of the transport equa- 
tions Eq. ([)]). The result is displayed in fig. ||, where we 
display the conductivity as a function of the electronic 
density for several temperatures. 




FIG. 5. The electrical conductivity, for the potential 
strength V z = 1.4 at different temperatures, calculated in 
EMA and with the Ziman formula. 

For each temperature the conductivity has a deep min- 
imum at n e ~ 1.8, corresponding to E m (see above). This 
minimum does not only follow from the short life time of 
the quasi-particles, but from large vertex corrections too. 
It becomes deeper and deeper as the temperature is low- 
ered. Note that the bundle of curves crosses at n e ~ 1.5 
and n e ~ 2.5, namely where the bundle of density of 
states crosses too. 

This is the root of the Mooij rule. Due to strong multi- 
ple scattering at the pseudo zone boundary, a pseudo-gap 
opens in the 1-particle spectrum, which becomes deeper 
and deeper as SRO grows. In turn it induces a min- 
imum in the conductivity which itself is getting deeper 
and deeper as SRO grows, namely as the temperature de- 
creases. It thus turns out that the resistivity decreases as 
the temperature increases, when the resistivity is large. 
Contrary to that, left from the minimum, the effect of the 



pseudogap results into the flattening of the effective dis- 
persion, leading to the maximum in the density of states, 
which becomes higher and higher as SRO grows, and so 
does the conductivity. It works however in the opposite 
direction right from the minimum, and as a result the 
resistivity can both increase or decrease as the tempera- 
ture increases, when the resistivity is small. 
The contributions to p (see eq. dTl])) are consisting of 
a contribution from the spectral function and a series of 
vertex corrections. It turns out that the latter are very 
importent, and can increase the contribution from the 
spectral function by as-jpauch as a factor 2. For a more 
detailed discussion seeEjlH3. 

In more detail the temperature dependence of the re- 
sistivity is displayed in fig. || for several values of the 
electronic density. In the regime of low electronic den- 
sity (n e < 1.7) where Ep < E m , the resistivity increases 
with temperature at a moderate rate. When the elec- 
tronic density n e approaches n m than the resistivity p 
decreases with increasing temperature T at a higher rate. 
It is maximum at n m . A further increase in the density 
leads to a decrease of this rate, for n e up to 2.4 where it 
changes its sign. 
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FIG. 6. The resistivity as function of temperature, for 
the potential strength V z = 1.4 and for different densities of 
electrons. 



Those data can be summarized by reporting the tem- 
perature coefficient of the resistivity a: 
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as a function of the resistivit 
est temperature leads to fig. 



Calculating a at the low- 
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FIG. 7. The Mooij-correlation for the potential strength 
V z = 1.4 at 300K, (o) corresponds to Ef < E m and (*) to 
Ef > E m - The minimum of a is correlated to Ef ~ -Em at 
n e = 1.78. 

The data scatter onto 2 lines which represent the 2 
regimes where the Fermi energy is either smaller or larger 
than E m . Not only the sign of a, but its amplitude as 
well, compare favorably with the exp erimentr 6 lr 3 . In con- 
trast to refJla we have calculated the temperature depen- 
dence of the resistivity for fixed density. There are 3 con- 
tributions to a: (i) the ^-dependence of the resistivity as 
such, (ii) the smearing of the Fermi surface together with 
a strong ^-dependence of cr(E), (iii) the change in the 
Fermi energy (at fixed density) following from a change 
in r\ and T. 

The first contribution to a is dominated by 1-particle 
effects and has been discussed above. This is the leading 
one. The second one results from the curvature of u{E) 
as one can infer it from the Sommerfeld expansion, and is 
contributing to a in much the same way as the first one; 
namely the curvature of <r(E) is positive in the vicinity 
of E m , negative below it, and very small above it. This is 
the second leading one. The third contribution, however, 
works in the opposite direction, and reduces the value of 
a at E m . 



C. Thermopower 



We now evaluate the thermopower according to Eq. 
( |l4| ) and we display it in fig. || where it is compared to 
the results obtained out of the Ziman formula. 

As compared to the results obtained out of the Ziman 
formula, we obtain much larger values, either positive if 
Ep < E m or negative if Ef > E m . Again the position of 
the Fermi energy with respect to E m plays the key role 
for assessing the sign of the thermopower. It changes its 
sign when the resistivity reaches its maximum value. 
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FIG. 8. The thermopower, for the potential strength 
V z = 1.4 and at different temperatures, calculated in EMA 
and with the Ziman formula. 

Thus phononic mechanisms are not required to ex- 
plain the experimentally observed large values of S, even 
though they are certainly going to play an important role 
as well, especially around the Debye temperature. An- 
other important contribution, namely the_phonon drag 
effect, is however small in metallic glassesEj, because the 
scattering of the phonons by the high degree of disorder 
keeps them essentially in local thermal equilibrium. 

We turn to the T-dependence of S and display it on 
ng. | 
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FIG. 9. The thermopower as function of temperature, for 
the potential strength V z = 1.4 and for different densities of 
electrons. 

Typically it is growing with the temperature when 
the Fermi energy is smaller than E m , whereas it is de- 
creasing when Ep > E m . If Ep is within fc^T of E m , 
there is a crossover region where dS/dT can have ei- 
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ther sign. Namely, if Ep < E m , we have dS/dT > 
if k B T < Ep — E m and dS/dT < if k B T > E F - E m . 
On the other hand, if Ep < E m , we have dS/dT < if 
k B T < Ep—E m and dS/dT > if k B T > E F -E m . Such 
a behavior has been observed in noble metals- polyvalent 
metals alloys such as (Au, Ag, Cu) - Sn. Indeed increas- 
ing the content of Sn (i.e. the density of charge carriers) 
changes the sign of the thermopower at low temperature 
from positive to negative while the resistivity decreases 
as review by HausslerE3. We finally note that the values 
of S are too large to fit the experiments. This follows 
from the relatively sharp E-dependence of the density of 
states, which is strongly reflected in the thermopower. 
The situation should improve by including multiple oc- 
cupancy corrections to the theory, but this is beyond the 
scope of this work. 



IV. SUMMARY AND ACKNOWLEDGMENT 

In summary we have solved the MT-EMA equations 
for a system of s-scatterers. The results lead us to a con- 
sistent interpretation of the anomalous transport proper- 
ties of amorphous metals. Especially we found that if the 
resistivity exceeds a certain threshold, than it decreases 
with increasing temperature, while otherwise the temper- 
ature coefficient of the resistivity can have either signs. 
This behavior follows from the opening of a pseudo-gap at 
the quasi zone boundary. The pseudo-gap is responsible 
for a large resistivity, as well as for a large thermopower 
(which can be positive or negative) and its non-trivial 
temperature dependence. As a function of the charge car- 
rier density, the low temperature thermopower changes 
its sign when the resistivity reaches its maximum. 
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V. APPENDIX 

In order to solve the set of equations Eq.(||) numeri- 
cally we use a set of transformations in order to make the 
task tractable. First we solve the equations for complex 
energies with Im(E) = 1 which turns out to be easy to 
solve. We then gradually reduce lm(E) down to 10~ 4 , 
a value that is typically smaller than — Im(E(fci?)) which 
has thus very little influence on the result. The equations 
will be solved by iteration. Even though we use a linear 
combination of the previous 3 iterations to improve the 
convergence, the combination being chosen so as to op- 
timize the difference between the input and the output, 
the task turn out to be very tedious and require as much 
as 10 3 iterations for each value of rj and Ep. We handle 



the on-shell singularity which enters the expression for 
G k (k) in the same way as described in ref namely we 
introduce: 



Gfc(«) 



G k ( K )-B k ( K ) 
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which obeys the equations: 
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(18) 



Nevertheless, even though the new equation is now 
smooth on-shell, it will develop a singular behavior at the 
renormalized quasi-particle peak through the T-matrix. 
The difference is that the location of the pole is not 
known beforehand, and lies in the complex plane. Thus 
some additional manipulations of the equations is re- 
quired. They are inspired by considering the integral: 



dx 



(19) 



The integrand is singular at z, z being complex. If Im(z) 
is large, the resulting integrand is smooth and such an 
integral can be evaluated as such. But when Im(z) is 
small, it is best to transform the integral Eq. ( ^9|) 



dx 



m - /(*) 



dx + f(z) 



1 



-dx (20) 



and evaluate the second term analytically. Obviously the 
resulting integrand is smooth. Applying this principle to 
Eq.(|j) allowed us to greatly improve the numerical sta- 
bility of the scheme. 

Another difficulty arises because the integrands in 
Eq.([|) decrease slowly for large momenta. In this case 
we have: 

1 - t(n, «)Ed + nt(n, k)B~ 1 (k)G ? (k) 2 + nt(K, k)G 9 (k) 
-B^~ 1 (k)(1 — t{n, n)(T,d + nGq(n))) — nt(K, k) 



a (k,q) 



2n 2 n 
Defining 



2n 2 n q 2 



H(k) 



-dq a (k,q) 



27T 2 7 



(21) 



(22) 
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we also obtain the large momenta contribution to as 



—nt{K, k) 

t(K,«)E d - 1 J k 



dk 

— {B k (n) + H{k)) (23) 



which is added to at each iteration, in the same way 

as H(k) is added to Gk{n). Both integrals entering Eq. 
(p2|, p3|) can be calculated once and for all. 
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